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Abstract 



We study the generation of density perturbations in the ekpyrotic scenario 
for the early universe, including gravitational backreaction. We expose inter- 
esting subtleties that apply to both inflationary and ekpyrotic models. Our 
analysis includes a detailed proposal of how the perturbations generated in a 
contracting phase may be matched across a 'bounce' to those in an expand- 
ing hot big bang phase. For the physical conditions relevant to the ekpyrotic 
scenario, we re-obtain our earlier result of a nearly scale-invariant spectrum 
of energy density perturbations. We find that the perturbation amplitude is 
typically small, as desired to match observation. 



1 



Typeset using REVTpX. 



We recently proposed a novel scenario for the early Universe in which the hot big bang is 
created by the collision between two M-theory branesS. The scenario assumes the Universe 
begins in an almost static, nearly BPS initial state consisting of empty, flat, parallel three- 
branes. In the effective 4d theory, the BPS state is homogeneous and has zero spatial 
curvature. Due to non-perturbative effects, however, a tiny force attracts the branes to one 
another. As the branes come together, quantum fluctuations create ripples in the brane 
surfaces that result in spatial variations in the time of collision. Consequently, some regions 
heat up and begin to cool before others, producing a spectrum of long wavelength density 
perturbations which can seed structure formation in the Universe. 

We estimated the perturbation spectrum using a 'time delay' formalism!, often used 
in simplified treatments of inflationary models. In that context, spatial variations in the 
time when inflation ends result in long wavelength density inhomogeneities. We applied 
the same formalism to variations in the time of collision in the ekpyrotic scenario. The 
equation for fluctuations in the scalar field <fi describing the inter-brane separation in the 
ekpyrotic model is almost identical to that describing fluctuations in the inflaton during 
slow-roll inflation. Consequently, a nearly scale-invariant spectrum of fluctuations is found. 
The result is remarkable because it shows that the Harrison-Zel'dovich spectrum can be 
obtained without inflation in a space-time which is very nearly static Minkowski space. 

The time delay formalism is a crude approximation, and only quantitatively accurate for 
a small class of inflationary potentials!. Nevertheless, it often gives a good estimate of the 
spectral index for the power spectrum of perturbations. One of the goals of this paper is to 
investigate whether the same statement is true for the ekpyrotic model. 

In the case of the ekpyrotic model, there is the major complication that the perturbations 
are produced when the effective 4d scale factor is contracting. In order to have a viable 
scenario, a mechanism must be found to reverse from contraction to expansion. This issue 
has been addressed in a recent paper we have written with N. Seiberg!, where we argue that 
such a 'bounce' may be allowed in the context of M-theory, where it corresponds to a collision 
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and rebound of the outer boundary branes. A matching rule linking the homogeneous 
background variables of the contracting phase to those describing the expanding phase was 
suggested there. Assuming this proposal is valid, what remains is to apply and extend those 
ideas to describe the evolution of perturbations through the moment of reversal. 

A stimulus for the present work was a paper by D. Lyth,! which calculated the growth of 
the perturbation variable ( (also commonly termed TV) representing the curvature perturba- 
tion of spatial slices which are comoving with the matter. Lyth correctly showed that £ was 
not amplified in the contracting phase of the ekpyrotic universe. He claimed this implied 
that when gravitational back-reaction was included, the spectrum of density perturbations 
became strongly scale-dependent with negligible power on large scales, making the ekpyrotic 
scenario incompatible with observations. His analysis employed a certain class of analyt- 
ically solvable models with exponential potentials, previously used to describe power law 
inflation! and simply extended to the situation of slow contraction relevant to the ekpyrotic 
scenario. We repeat his analysis here, but also compute the perturbation in the Newtonian 
potential $. We show that in the contracting phase gravitational back-reaction actually en- 
hances rather than suppresses long wavelength fluctuations, but these fluctuations show up 
purely in $ and not in (. Therefore gravitational back- reaction does not spoil the ekpyrotic 
mechanism, at least in the contracting phase. 

The remaining issue regards the appropriate matching condition for tracking the per- 
turbations across the bounce and into the expanding hot big bang. Consistent with the 
arguments of Ref. 4, we seek to identify variables which are non-singular at the bounce, 
both for the background and perturbation variables. We then match the amplitudes of the 
two linearly independent solutions for the perturbation variables across the bounce. With 
our prescription, we find that the long wavelength perturbations developed in the contracting 
phase do indeed survive to the expanding phase, provided there is a change in the equation 
of state at the bounce, such as occurs if a sub-dominant component of radiation is produced 
there. Our final expression for the density perturbation spectrum agrees well with the more 
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naive time delay estimate. 

After communicating a preliminary version of this paper to Lyth and R. Brandenberger, 
Lyth prepared a second manuscript^ proposing that contraction be matched to expansion 
on a time-slice of fixed energy density^. With this procedure, he argues that the curvature 
perturbation ( is conserved across the bounce. But since ( does not acquire a scale invariant 
spectrum in the contracting phase, he argues that ( will not have such a spectrum in the 
expanding phase, where it represents the amplitude of growing mode adiabatic perturba- 
tions. Lyth's conclusion is that any growing mode density perturbations developed in the 
contracting phase match perfectly onto pure decaying mode perturbations in the expanding 
phase. Brandenberger and Finelli,! and Hwang,§ have recently produced preprints repeating 
this argument. In a note added, at the end of this paper, we explain why we do not believe 
these conclusions are valid for the ekpyrotic scenarios proposed in Refs. 1 and 4. 

Let us outline our approach to the matching problem. We want to evolve background 
and perturbation variables according to the appropriate field equations, all the way to to 
zero scale factor in the four dimensional effective theory. We identify a complete set of 
variables which are non-singular at the bounce, and match those non-singular variables 
across it. This prescription automatically excludes the variables $ and (, both of which 
diverge. More generally, geometrical quantities such as the synchronous gauge comoving 
metric perturbation hij also diverge. Indeed the meaning of the three-geometry is unclear 
at zero scale factor. 

Instead our approach is essentially algebraic rather than geometrical. We focus on gauge 
invariant perturbation variables which are consistently small at all times and match these 
across the bounce at t = 0. We argue that this matching would give consistent results for an 
infinite class of perturbation variables so defined. Our prescription can only be fully justified 
by a satisfactory microscopic description of the relevant degrees of freedom. Nevertheless, 
if string theory shows that the scale factor can truly pass through zero and bounce, then 
tracking perturbative gauge- invariant degrees of freedom which remain small and finite seems 
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likely to be the right approach to matching fluctuations across the bounce. 

Our analysis shall be performed entirely within the context of four dimensional effective 
field theory. This does not capture all the low energy degrees of freedom relevant to the 
five (or indeed eleven) dimensional brane- world. As in Ref. 1, we shall assume those other 
degrees of freedom are frozen, or at least so slowly varying that their inclusion would not 
substantially alter the result. We shall focus here on single moduli field which determines 
the outer brane separation in a brane-world Universe. In Ref. 1 we considered a model 
in which the perturbations are produced by the collision of a bulk brane with one of the 
boundary branes. As we pointed out there, if the 4d effective theory is valid, the scale factor 
in that theory must continue to contract until the outer branes collide and bounce. In Ref. 
4 we suggested a simplified model in which there is only one collision between the boundary 
branes and no bulk brane is needed. The perturbations are produced as the outer branes 
approach one another. For simplicity, here we shall restrict ourselves to this two-brane 
scenario, in which the same scalar field (f> is responsible both for the development of the 
perturbations, and for describing the the outer-brane collision and bounce. Generalizations 
to bulk-boundary collisions, as in the original scenario of Ref. 1, are a straightforward 
extension and will be briefly mentioned when appropriate. 

The outline of this paper is as follows. In Section I we discuss the properties of the inter- 
brane potential relevant to the ekpyrotic scenario, an issue referred to again in the note 
added at the end of this paper. We then review the application of the time delay formalism 
to the ekpyrotic scenario, as described in Ref. 1, showing that a scale invariant spectrum of 
fluctuations is naturally predicted. In Section II, we show that including gravitational back- 
reaction has, in Newtonian gauge, only negligible effects on the fluctuations acquired by the 
scalar field. Section III is devoted to a discussion of the role of the curvature perturbation £ 
(or TV) conventionally used in the analysis of inflationary models, which is in fact insensitive 
to the growing mode perturbation in the contracting phase of the ekpyrotic model. This 
is further elaborated in Section IV where we show that ( is canonically conjugate to the 
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direction of amplification, which is proportional to the variable $. As a consequence of 
Liouville's theorem, ( is 'squeezed' as the perturbations develop during the collapsing phase 
of the ekpyrotic scenario. 

Section V addresses the key moment of reversal, which is when, in the four dimensional 
effective description, the Universe 'bounces' as the scale factor a hits zero. Since the four 
dimensional spacetime geometry is singular there, a four dimensional geometrical description 
using the Einstein frame metric is inappropriate. However, our approach is to identify a 
complete set of dynamical variables which remain finite at the bounce and to match them at 
the transition from contraction to expansion. We ignore divergent quantities as unphysical, 
attributing their bad behavior to a singular choice of field and metric variables. We used this 
approach in Ref. f| to describe the homogeneous background evolution during the bounce. 
Here we shall show that there is also a set of gauge invariant linear perturbation variables 
which are well-behaved at the bounce. When the theory is formulated in these variables, 
there is a well defined matching prescription, even in the four dimensional effective field 
theory, which seems to yield physically sensible results. 

In order for the well-behaved variables to exist, we emphasize that it is important to 
have two conditions: (1) there must be a free scalar field (the modulus in our case) whose 
kinetic energy diverges at the bounce and (2) the scalar potential V{4>) must be sufficiently 
non-singular at the bounce. The simplest case is where V(4>) approaches zero at the bounce, 
as suggested by the mapping from M-theory to weakly coupled string theoryi. If these two 
conditions are fulfilled, then it is possible to follow the perturbations and match at a = 0. 

Our main finding, in Sections V and VI, is that as long as radiation is produced or fields 
are excited at the outer-brane collision so that there is a jump in the first or second time 
derivative of the equation of state parameter w = P/p, the scale invariant perturbation 
spectrum developed during the contracting phase propagates through the bounce and into 
the final expanding Universe. Intriguingly, the final density perturbation amplitude derived 
in Section VI is naturally suppressed by a small numerical coefficient, by quantities which 
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are automatically small in Horava-Witten theory, and by factors involving the efficiency with 
which brane kinetic energy is converted into radiation. We comment on how this suppression 
may naturally explain the small perturbation amplitude ~ 10~ 5 we see via observations of 
the cosmic microwave sky in the universe today. 

I. TIME DELAY NEGLECTING GRAVITY PRE-COLLISION 

We first derive the ekpyrotic spectrum in a very naive 'time delay' approach which totally 
neglects cosmological expansion, gravity and all moduli other than 0. We also ignore the 
crucial element of reversal from contraction to expansion, which will be a critical aspect of 
the discussion in Sections V and VI. As we shall see, despite the fact that the time delay 
argument ignores these very important features, it nevertheless comes close to matching our 
final answer. Thus, as is the case in inflation, the time delay argument turns out to be a 
convenient heuristic even though it is neither physically rigorous nor numerically accurate. 

We assume that at large positive the effective potential governing the evolution of 
takes the form 

V = -V e- C <t>, (1) 

where Vq and c are positive constants. The ekpyrotic scenario starts at large positive 0, in 
a state of nearly zero energy, and rolls towards negative values. As becomes increasingly 
negative and the branes come close together, V(<f>) must turn upwards and approach zero. As 
discussed in Ref. |], for example, this is the behavior to be expected in M-theory because the 
string coupling constant vanishes as the outer branes collide. (A similar constraint applies 
in the bulk brane collision, as discussed in Ref. 1.) 

Generally, in brane world models, the separation of the two boundary branes can be 
described by a canonically normalized, minimally coupled scalar field (the 'radion'). As 
the brane separation d goes to zero, one can ignore the effect of the bulk 'warp factor' and 
one obtains the Kaluza-Klein result in which d ~ e \/V^'t>/ M Pi q as __>. _ oo. (Here and 
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below, Mpi = (87rG) _ 2 is the four dimensional reduced Planck mass.) In standard Kaluza- 
Klein theory, the same dependence holds at large 0, so the range of is — oo < < oo. 
But the presence of a bulk warp factor alters this. In many models, when the outer-brane 
distance d tends to infinity, tends to a finite value, which may be taken to be zero. For 
example, for a pair of positive and negative tension boundary branes with a bulk Anti-de 
Sitter space, the inter-brane distance is d = L ln(coth(— <p/^/6Mpi)) where L is the AdS 
radiusi. A suitable ekpyrotic potential would then behave as ~ — (— <p) N ~ — e^~ Nd ^ L \ at 
small 0, and we would be interested in starting the system in a state such that — > at 
large negative times. In models with a bulk brane, such as those considered in Ref. 1, the 
location of the bulk brane would again be described by a scalar field, but this time its range 
would be finite. The exponential model (|1|) is useful, since it is mathematically tractable. 
However, one should remember, in general it should only be expected to apply over some 
finite range of 0. 

We now review how the time delay formalism can be applied to our example. The 
classical solution is given by solving = — y— 2V. For V((f>) given in (J]]), we obtain 



where the time t is large and negative at large 0. We follow the evolution to some finite, 
small negative t, at which point the modes of interest are 'frozen in'. For small fluctuations, 
expanding in plane waves 50 = X^ eJfc x <50fc(^O; one finds 



static, this is the same equation as that governing a massless field x in a de Sitter background 
(with x — 0/ a ; with a the scale factor and t conformal time). As in that case, starting from 
an initial Minkowski vacuum one generates a long wavelength spectrum of scale invariant 
fluctuations. 



k 1 , so the amplitude of normalized modes with wavenumber k is ~ k 1 ^ 2 . Consider a mode 





<% = -/c 2 <% - V M 5^ k = -kH(j> % + -% (3) 
(using the solution Eq. @ above). Curiously, even though the background here is nearly 



We assume that the quantum fluctuation 50£ starts in the Minkowski vacuum, (|50^| 2 ) oc 
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where — kU 3> 1, where U is the initial time. The mode oscillates at fixed amplitude 
k~ l/2 until t increases to t\ ~ —k l , when the destabilizing 2/t 2 term in (|3|) begins to 
dominate over the stabilizing k 2 term. One then has S(j)^ ~ (2/t 2 )<5</>£, with growing solution 
5(p£ oc (— t)" 1 . Hence, the final amplitude for the fluctuation, as t approaches zero, is 
<50£ ~ k~ l ' 2 {— t) ~ /c~ 3 / 2 (— t)~ x . The precise result for the long wavelength spectrum 
is 

<i%r> = 2^ -kt«i, (4) 

which is scale-invariant as claimed. 

A naive estimate of the final density perturbation amplitude runs as follows. The scalar 
field fluctuations give rise to a time delay in the moment of collision. Neglecting gravitational 
effects, the time delay between collisions relative to the mean collision time is — <50/0o- The 
net density perturbation after collision is then just y ~ —AH5t where H is the Hubble 
parameter at the time of collision. Despite the deficiencies of this time delay approach, 
this answer is quite close to the result we shall eventually derive, including gravitational 
back-reaction and matching across the bounce, which we give in equation ( |B1D below. 

The first step in improving the above treatment is to include gravitational backreaction 
in the initial, contracting phase. For this purpose, we focus on the gravitational backreaction 
corrections to Eq. (H). 

II. INCLUDING GRAVITATIONAL BACKREACTION 

For exponential potentials, there exist analytic scaling solutions to the Friedmann- 
Robert son- Walker equations which at large negative t approximate the assumed initial con- 
ditions in the ekpyrotic setup, and which allow an analytic treatment^. For the exponential 
potential considered as an approximation over some range of in the ekpyrotic scenario, 
one has the background solution 

a{t) = (-ty, Mt) = -H-Mt), (5) 

c 
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where a(t) is the scale factor, t is proper FRW time, taking negative values, and 

P ~ - M>(l-3p)' W 

and Mpi = (8nG)~^ is the reduced Planck mass. Some useful quantities are: 



where dot denotes d/di. In ekpyrotic models, we are interested in potentials in which the 
evolution is very slow or, equivalently, for which p « 1. (Parenthetically, note that the same 
solution, for p > 1, describes a universe undergoing power law inflation). The conformal 
time is 

fo dt t 



(8) 



It a(t) (1 — p) a(l — p) ' 
running from — oo to as t does, although as mentioned we shall have to alter the form of 
the potential as t and r approach zero, to avoid its diverging to minus infinity. Derivatives 
with respect to t shall, as above, be denoted by dots and henceforth, derivatives with respect 
to r shall be denoted by primes. 

Mathematically, Newtonian gauge is convenient for studying scalar field perturbations 
since the linearized equations reduce to a single second order differential equation for the 



Newtonian potential $ (for a derivation see e.g. Ref. [10]) 



$" + 2 ^$' + k 2 $ + 2<j)' Q (^j $ = 0, (9) 

where A = (0o) _1 , an d TC = a' /a. Here, and below, the wavenumber dependence of all 
perturbation quantities is not shown explicitly. For $ one should read <3>£ and so on. 
We now eliminate the first derivative term by setting $ = u/A, obtaining 

A" ( H \ ' 

u" - —u + k 2 u + 20o U- J u = 0. (10) 

A V0o/ 
Substituting the above scaling solution, we obtain 

u " = ~ k2u + (I^}V^ ( U ) 
10 



which is a form of Bessel's equation, with corresponding order v = / (1— p). Note that 

our notation for the variable u, and for the variable v introduced later, matches Mukhanov's 
original notation^. 

The initial conditions are that the scalar field fluctuations should be in the Minkowski 
vacuum state as r — > — oo: in Newtonian gauge, two constraints determine $ and $ in terms 
of 86 and 86: 



2M|, 



2M 2 Pl 



^+r- + 6 d t [^) }9 



k 2 , „ H 



a 



-iki 



(12) 

At large — kt, for 86 ~ e - %kr j (a\/lk) , which is the incoming Minkowski vacuum, one finds 

$ „ ^ e -ifcr (13) 

2MpikH 

which vanishes at early times, consistent with our ekpyrotic initial condition that the space 
is asymptotically Minkowski in the past. 

Neglecting an irrelevant phase factor, the solution for u is then 

u=(-kr)Kl^Hl(-kT) ^ , (14) 

1 V 2 ; (2k)lM Pl 1 ; 

where if* is a Hankel function and v was given above. We follow this solution forward to 
conformal times at which — kr << 1, when the modes become frozen in. Converting back to 
$, and using the small argument expansion for the Hankel function, again neglecting phase 
factors we find 

$ ~2^^ (1+I,)( " rr2 ' (1 + ° (P)) - (15) 

Since we are interested in p << 1, we henceforth neglect corrections of order p to the 
numerical coefficient. But we keep the p dependence in the scaling with wavelength and 
time for comparison with later results. 

Now we can convert back to the scalar field 86, using the formulae ( ^ ) given above, and 
obtain: 
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9 M 2 r 

^ = ZMpi ^ + ^ 2 - 1 /2 A ;-a+-)(-t)-(l+f). (16) 

where we used $ oc (-r)- 2 ^ oc (-t)-( 1+p \ The resulting power spectrum of fluctuations is: 



where {\5(j)£\ 2 )^ is the result obtained in Eq. (d) ignoring gravitational backreaction. Re- 
calling that the regime of interest for the ekpyrotic scenario is small p, we see that we obtain 
the same answer as in Eq. (H) up to small corrections. The backreaction produces a slight 
reddening of the power spectrum, presumably because long wavelength perturbations were 
generated sooner and had longer to self-gravitate. However, since in the ekpyrotic scenario 
we have p « 1, these corrections to the spectral index are actually smaller than corrections 
arising from the nontrivial kinetic term for <f), which were computed in Ref. 1. The result 
shows that backreaction and metric fluctuations have an insignificant effect on the density 
perturbations developed during the initial, contracting phase of the ekpyrotic universe. 

III. THE CURVATURE PERTURBATION IN THE CONTRACTING PHASE 

Rather than solve for the evolution of the Newtonian potential in Eq. @, a common 
procedure used for inflationary models is to track the curvature perturbation on spatial slices 
which are comoving with the matter. As we shall discuss, this variable is in fact insensitive 
to the growing mode density perturbation in the contracting phase. 

Following the notation of Mukhanov et a/., Ill we denote the curvature perturbation on 
comoving slices by (, defined by 

C = 3^T^ + * = ^G^ + *' (18) 
where 7i = a' /a and w is the ratio of pressure to density in the background universe. The 
variable £ was introduced by Bardeen0 in his classic paper, in which it was termed <p m . It 
was employed in the context of inflation by Bardeen, Steinhardt, and Turner^ and its use 
later elaborated upon by many authorsilllS. Mukhanov et alM re-express ( as ( = v/z 
(for spatially flat hypersurfaces) and thereby derive the 'w-equation': 
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z" 



v " = -k 2 v + — v (19) 

z 

where z = a(p' /H. For the power law solution given in equation (|7|), this reads: 

Since the last term is negative, it follows that there is no classical instability in the variable 
v, a point emphasized by Ref. 3. (Note that Ref. 3 relabels Mukhanov's v variable as u; here 
we have kept to the original notation.) 

However, this does not at all imply the absence of a growing mode density perturbation: 
the variable $ does exhibit an instability, as we have discussed above. Given a solution v, 
the potential $ can be obtained by solving 

k 2 <$> = -AnG^v/z)'. (21) 

The general solution to the f-equation to order k 2 is: 

v = C lZ (l-e T % f ' z 2 dr") + C 2 z( f %- k 2 [ % f ' z 2 dr") , (22) 

V J —oo Z J — oo / V J— oo Z J —oo Z J —oo I 

where C\ and C 2 are arbitrary constants. Substituting into Eq. (Ell), we obtain 

$ = C x % I" z 2 dr' + C 2 <!> (23) 



where $o = of /a 3 , and Co and C 2 are constants. For the power law solutions studied here, 
z oc a. For an expanding universe, the first term dominates the second as t oo. However, 
for a contracting universe, the second term is the proper growing mode as t — > 0. The 
solution is the same as that found in Eq. (|i5|), from which we derived the scale-invariant 
spectrum. 

Various subtleties are worthy of further comment. First note that $o = of /a 3 is an exact 
solution of the perturbation equation (|9|) for k = 0. In fact, for k — 0, the perturbation $o 
just represents a coordinate transformation of the the time t, and is therefore unphysical. 
However, for nonzero k the solutions tending towards $o are not gauge modes and therefore 
represent a real gravitational instability. 
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Second, note that ( is insensitive to any component of <3> tending towards $o> as may be 
seen from the expression 



C 3a 2 (l + w) (a'/a 3 J ' (24) 
For inflationary models, tracking (, as is done in many analyses, is useful because the 
components that are projected out in this formula are harmless decaying modes. However, 
for the ekpyrotic model, the growing mode tends to $o- O ne can in fact recover the the 
growing mode solution from ( or v, but special care has to be taken to do so, by keeping 
the next order terms in k 2 as we have done above. 

IV. CANONICAL CONJUGATES AND 'SQUEEZING' 

In this section we show that the variables u in Eq. (|T0| ) and v in Eq. at the heart 
of the analysis of the previous sections are in fact canonically conjugate. The importance of 
this is that in the (p, q) phase space plane, the trajectories of the system are focused towards 
certain lines - for example the line p = q in the case of the upside-down harmonic oscillator. 
These lines are then the classically amplified directions in phase space, and from Liouville's 
theorem the phase space density is necessarily squeezed in the orthogonal directions. While 
the calculation of Reference 5 is indeed correct, it does not by itself indicate the absence 
of a growing mode density perturbation, because the curvature perturbation which was 
computed is actually orthogonal to the direction in which fluctuations are amplified in the 
ekpyrotic setup. (We note that References [7| and |8| agree with this conclusion). 

Our starting point is the action for gravity plus a scalar field in canonical (first order) 



form. The relevant formulae may be found for example in Ref. [15]. The Newtonian potential 
$ is written in terms of u, and ( is written in terms of v as above. One finds the quadratic 
action for perturbations reduces to 

S = 2 -Mk ( dr (vF' + -vF - -F 2 - —v 2 ) , (25) 



p J \ t 2 2 

where (3 = p/(l — p), F = k 2 u, and we have written the Laplacian V 2 as — k 2 . The sum 

over all Fourier modes is implicit. 
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The action ( p5|) is in canonical form J pg' — H(p, q), with if the Hamiltonian. Therefore, 
up to normalization, v and F are canonically conjugate. The equations of motion are 



y k 2 v = r- p (FT P )', (26) 



which are just the relations above: equation ([H]), which expresses $ in terms of ( and 
and equation (|24f) , which expresses £ in terms of $ and These imply the equations of 



motion (|TT|) and (|T^) used in the previous sections. 

In the inflationary case, which corresponds to p > 1 in the above analysis, both the 
variable v and its canonical conjugate momentum F (or u) obey equations exhibiting a 
classical instability: both are amplified as r approaches zero. The system starts near v = 
u = 0, with only small quantum fluctuations. As time proceeds the system evolves away 
from the origin along a certain line, u cx v in phase space. By Liouville's theorem, the phase 
space density is 'squeezed' in the direction orthogonal to this line. But in the ekpyrotic case, 
p « 1, the v (or £) variable is instead stable and classical trajectories instead run out along 
the u (or $) axis. The v (or Q direction must in consequence be 'squeezed'. Thus there is 
no classical growth in (, consistent with the finding of Reference |5|. 

V. REVERSAL AND A NONSINGULAR MATCHING CONDITION 

Having established that there is indeed a growing mode scale-invariant spectrum of den- 
sity perturbations in the collapsing phase, we now turn to the more challenging question of 
whether these perturbations survive the reversal to expansion. We have recently addressed 
the issue of reversal in a paper with N. Seiberg,§ and we shall employ and extend the analysis 
of that paper here. 

As we have discussed, the variable ( is not amplified in the collapsing phase, since it is 
insensitive to the growing mode perturbation. However, ( does yield the amplitude of the 
growing mode perturbation in the expanding phase. Therefore, if £ was in fact the correct 
variable to match at r = 0, the growing mode perturbation which developed in the collapsing 
phase would just match onto a pure decaying mode perturbation in the expanding phase. 
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This has indeed been suggested by Lyth in a follow-up paper following our communication of 
the above calculations to him0. The new argument is essentially that perturbation evolution 
may reverse at collision. That is, $ began negligibly small, grew to be very large at the 
bounce but then, after reversal to expansion, $ shrank precisely back to zero. 

In this section, we wish to explain why matching ( across the bounce is not appropriate 
in the scenarios we consider. Instead, we claim, the scale invariant growing mode spectrum 
generated in the contracting phase matches onto a linear combination of growing and de- 
caying modes at reversal. Consequently, the density perturbations survive the reversal of 
the scale factor to seed structure formation in the hot big bang. 

In order to analyze matching at the bounce, gauge invariant variables must be identified 
which are well-behaved as a and r approach zero. We have already seen that the growing 
mode solution for $ is proportional to a' /a 3 , which is divergent as r tends to zero. Likewise, 
as discussed below, in the situation of interest, ( is logarithmically divergent. Since linear 
perturbation theory can only be valid when the perturbation variables are small, it follows 
that ( is not a good matching variable. Nevertheless, with careful treatment to exclude 
the logarithmic divergence, ( shall be very useful in the analysis after r = 0, just because 
it's 'long wavelength' component is nearly constant in the expanding phase, and gives the 
amplitude of the growing mode linear density perturbation. 

As mentioned in the introduction, the potential V(4>) in (HD cannot hold as <fi — > — oo. In 
the main example we treat here, we instead assume V(<fr) bends upward so that V(<f)) -»0as 
— > — oo. In this case, just before collision the system is completely described as a massless 
scalar field coupled to gravity, in a collapsing Universe. As discussed in Ref. 4, one can in 
this case define regular variables, in which the evolution through a = becomes well defined. 

For the perturbation analysis, it is very important that, as a — > 0, the Universe is dom- 
inated by scalar field kinetic energy. In this situation, a good perturbation variable is the 
fractional energy density perturbation on spatial slices which are comoving with the mat- 
ter, termed e m by Bardeen in his discussion of gauge invariant cosmological perturbation 
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theory .0 This variable (defined in his equation 3.13) is a linear combination of the energy 
density perturbation 5T®, and the velocity perturbation T?, which is invariant under lin- 
earized coordinate transformations. It equals the energy density perturbation in any gauge 
in which the matter worldlines are orthogonal to the r =constant hypersurfaces, and as 
Bardeen emphasized it is the natural choice of perturbation amplitude from the point of 
view of the matter. 

The equations of motion for the matter, = 0, lead to the following equations of 
motion, 

(pa 3 e m )' = -(p + p)a 3 kv^ 

vf + Kv® = M> + y-^fe + w V ), (27) 

(1 + w) 

where is the gauge-invariant velocity perturbation defined by Bardeen. Here we intro- 
duce w = p/p parameterizing the equation of state for pressure p and energy density p, 
c 2 s = dp I dp being the sound speed, and rj being the 'entropy perturbation', defined as the 
difference between the pressure perturbation and that expected from the density perturba- 
tion and the background pressure-density relation. In most of the analysis below we shall 
assume rj is zero, which is to say that the density perturbations are 'adiabatic'. For the 
types of matter we consider (perfect fluids and scalar fields), the anisotropic stress is zero, 
and Bardeen's potentials are related to our potential $ by $ = $^4 = — With these 
simplifications, and for the flat universe we consider, equations (4.5) and (4.8) of Bardeen's 
paperS reduce to ( [27]) above. 



Remarkably, scalar field kinetic energy domination implies that e m is actually finite at 
t = 0. This can be seen from the Einstein constraint equation, which reads 

V 2 $ = A7iGpa 2 e mj (28) 

where p is the total background energy density. For scalar field kinetic domination, p oc a -6 
as a tends to zero. But the growing mode perturbation $ oc a' /a 3 oc a -4 and, from Eq. 



), e m is finite as a tends to zero. 
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Having identified a non-singular matching variable, we still need to decide which time 
slice to match on. In the present situation, this is unambiguously defined as time slice 
where tends to — oo, or, re-phrased in terms of the brane separation, the time when 
d oc e ^^^ Mpl — > 0. If we assume that no radiation is present in the contracting phase, so 
that the only energy density present is that in 0, then the surfaces of constant are also 
the comoving surfaces (because the perturbation to the momentum is proportional to S(f>), 
and e m is the fractional energy density perturbation evaluated upon these surfaces. 

At the bounce, we assume that there is a change in the internal state of the branes so 
that V{4>) switches off and the boundary branes no longer attract. After the bounce, is 
therefore a massless, free field. We also assume some radiation is produced on the branes. 
Again the matching surface is defined in terms of 0, but with radiation present it is not 
obvious that the surfaces of constant are still comoving with the matter. However this is 
indeed the case for adiabatic perturbations as the following calculation reveals. 

We perform the calculation in conformal Newtonian gauge. The equation governing the 
evolution of the radiation perturbation is given in this gauge by 

(5 r -4$)' = hv r , (29) 

where 5 r is the fractional energy density perturbation in the radiation and v r the velocity 
perturbation. We now make the assumption of adiabaticity, equivalent to the statement that 
the ratio of radiation to scalar field energy density is spatially uniform being determined 
locally by the microphysics of brane collision. The condition for the 'entropy perturbation' 
rj, defined above, to be zero at early times is that 5 r = |<^,, where 5^ is the fractional density 
perturbation in the scalar field. In Newtonian gauge we have 

*=2(£-.). (30) 

The final ingredient in calculating v r is the equation for perturbations in the massless scalar 
field, 

5(f)" + 2H5<f)' = 40(,$' - k 2 5<f). (31) 
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^From the adiabaticity condition and equations (p9j), (|30|) and (pTj) , as well as the background 
equation 0q + 27i<// = 0, one determines 



This is the velocity of the radiation in Newtonian gauge. We wish to transform to the time- 
orthogonal gauge in which 5<p is zero. The transformation required is a shift in conformal 
time by T = —S(p/(j)' accompanied by a shift in spatial coordinates with a potential L 
such that L' = —kT, to ensure the absence of time-space components of the metric (in the 
notation of Section III in Bardeen's paper0). The result of these transformations is that 
in the 5(j) = gauge, the velocity potential for the radiation equals v r + L' , which is zero. 
Therefore, for adiabatic initial conditions, the radiation fluid is actually comoving with the 
massless scalar field. Therefore the comoving slices correspond to constant scalar field slices, 
both before and after the bounce. 

So we may therefore focus on the evolution of the energy density perturbation on comov- 
ing slices, e m , and attempt to match it across the bounce at r = 0. The equation of motion 
for e m is straightforwardly derived from ( |2"T| ) by eliminating v^°\ to obtain 




(32) 




(33) 



where 



f(r) 



H(l + 3c 2 s - 6w) 




(34) 



As an aside, we now infer the behavior of the variable ( at small r. One has 



c 




(35) 



Employing the background relations 





(36) 
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one obtains 

Hc 2 e 
1 + w 

Since (' is down by k 2 relative to (, the rate of change of £ is typically much smaller than 
7i( for the low k modes of interest. So one might imagine ( would be a good matching 
variable. However, as noted above, e m is finite as r tends to zero. From Eq. fl3T|) and 
the fact that Ji oc r _1 for small r, one observes that ( diverges logarithmically at small r. 
Hence, as stated above, in linear perturbation theory ( is not a good variable for establishing 
a matching condition at r = 0. Nevertheless, ( does have some utility. If one separates out 
the short wavelength piece which is divergent, the remaining long wavelength component 
of ( is finite and nearly constant. This finite piece is then very useful, since it gives the 
amplitude of the final growing mode linear density perturbation in the expanding Universe. 

We now turn to an analysis of Eq. (j33l). The first thing to note is that in the setup 
considered here, 7i, w and (? s all have expansions in terms of simple powers of r, for either 
r < or t > 0: 

H = ^r" 1 + h + h 1 r + 

W = 1 + W\T + W2T 2 + . . . , 

c 2 s = 1 + cit + c 2 r 2 + .... (38) 

This follows from computing 

<H la2 ' 

_ \<P 12 - a 2 V _ 3M 2 Pl Q 2 - a 2 U 
W ~ \<P 12 + a 2 V ~ ZM 2 Pl Q 2 + a 2 U' 
2,2 Vtpa 5 _ 2 a l U fi + a Q U il 

° s + 3a'0'a 2 + 3 HM Pl Q ' 1 ' 

where a 2 = \{oq — af), U = Va 4 and Q = ^(aoa[ — aia' Q ). As discussed in Ref. 4, the variables 
ao = 2a cosh(0/ v6Mp;), and a\ = —2asmh((f)/\/QMpi) are both finite at r = 0. These 
variables were introduced in Ref. [TJ, with motivation from the five dimensional geometry. In 
the AdS case, a® and a\ are the scale factors on the positive and negative tension boundary 
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branes respectively. Although the four dimensional Einstein frame scale factor a tends to 

zero at an outer-boundary collision, <f> tends to minus infinity in just such a manner as to 

leave a® and a\ finite at collision. 

The physics of bounce and reversal are discussed in Ref. |], the results of which are 

assumed here. The key point for our analysis is that the equations of motion for a and a 1; 

3 , r o „ dU 3 , , 9 „ dU , s 

2 M " a ° = ft£ 2 M " a " = -8ST- (40) 



are actually regular at r = 0. It follows that all the variables in Eq. (p9[) have simple power 
series expansions around r = 0, as shown in Eq. ([38]), for either positive or negative r. 

It will be important below to include a component of radiation after the bounce. As 
discussed in Ref. f|, the density and temperature of radiation on the branes is actually finite 
at a = 0, as long as the radiation couples to the scale factors a or a\, which are both finite 
(and equal) as <p —oo. If, for example, U is zero after the collision but there is a density 
p r of radiation produced on the brane with scale factor a (in the simplest models, we have 
a = ao or a% for the positive or negative tension branes respectively). Then the expression 
for w in Eq. (|39| ) is replaced by 

3M^Q 2 + l(p r aV 
W 3M^Q 2 + (p r a 4 )a 2 ' 1 ' 

Before proceeding to analyze Eq. (|33|), it is important to recognize that the expansion 

coefficients in Eq. ( |38D are not independent, but are related by Eq. (|36|) . These imply, for 

example, that ho = — j^Wi and c\ = The coefficient functions in Eq. (|33]) have the 

following expansions 

/(r) = -- + /o + /iT + ..., g(r) = g^r' 1 + g + 9l r + . . . , (42) 
r 

with /o = g_i = and g = ^wf + k 2 . Now the general solution of the 

perturbation equation (Eq. (|33"D) is a sum of the two linearly independent solutions 

e m = e D(r) + e 2 E(r), (43) 

where eo and e 2 are arbitrary constants, and 
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D(t) = l + diT + d 2 r 2 \n\r\ + d 3 T 3 ln|r| + d 3 r 3 + 

E(t) = r 2 + e 3 r 3 + . . . , (44) 

The equation of motion determines all the coefficients d n , n = 1, 3, 4, . . ., d n , n = 2, 3, . . . and 
e n , n — 3, 4, — One finds d\ = g_\ = jWi, d 2 = — |(di(/o + 9-i) + fl'o) = ~~ \k 2 , and so on. 
Since the coefficients w±, w 2 , ■ ■ ■ change across r = 0, the two series expansions are different 
for r < and r > 0. We denote the coefficients in (f43"[), for r < and r > respectively, as 
e (0~), e 2 (0~), and e (0 + ), e 2 (0 + ). The matching rule we seek should determine the latter 
two constants in terms of the former. 

It seems clear that we should match the amplitude of the finite perturbation variable 
e m across r = 0, and this fixes eo(0 + ) = e (CT). But how should we determine €2(0 + )? 
The simplest prescription is just to set 62(0 + ) = 62(C) - ). This amounts to matching the 
amplitude of the linearly independent solution which vanishes at r = 0, as well as that 
which is finite at r = 0. This prescription is invariant under redefining the independent 
solutions, e.g. by adding an arbitrary amount of the solution E(r) to D(t). Matching any 
other non-singular perturbation variable, defined to be an arbitrary linear combination of 
e m and e' m with coefficients which are non-singular background variables (defined to possess 
power series expansions in r, as above) will, with the same prescription of matching the 
amplitudes of both linearly independent solutions, also yield precisely the same result. 

This prescription is simple, but it is certainly not unique, and we emphasize that a 
proper understanding of the correct matching condition must ultimately rest on a better 
understanding of the singularity, either directly from string or M theory or from a well 
defined regularization and renormalization procedure. 

Note that any such matching rule applied at r = cannot match both e m and e' m , as 
would be appropriate at a regular point of the differential equation. This is because e m 
and e' m cannot be independently specified at r = 0, just because r = is a singular point 
of the differential equation. The value of e 2 is fixed by the limit as r tends to zero of 
\{.£m(j) ~ £qD{t))" ', so one needs to know e" m as r tends to zero, in order to determine e 2 . 
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Recall that since ( diverges logarithmically as r — > 0, it is not a good matching variable. 
However, this divergence only affects the short wavelength part of (, and the long wavelength 
part is still very useful since it yields the amplitude of the growing mode perturbations in the 
expanding phase. As we now see, the above prescription for matching e m actually implies 
that the long wavelength part of ( has a jump across r = 0. Using Eq. (^), we re-express 
Eq. (|D as 

c = - (TTSjF ( 4 + 5 H(1 " t " )e "') ' (45) 

Substituting the above expansions, one finds the leading order behavior 

C~f(ln|r|-^) + ^6 2 -A eoW ( 2 ) ); (46) 

where it/ 2 ** = u>2 + § w i; phis terms which vanish as r tends to zero. The first term is 
logarithmically divergent at r = 0. However, since it is down by a factor of k 2 , it rapidly 
becomes irrelevant as r increases away from zero. The second term, the long wavelength 
piece ( lw , is the quantity we are actually interested in. This constant, long wavelength piece 
is accurately conserved after the bounce as long as the matter evolution remains adiabatic, 
and yields the amplitude of the growing mode adiabatic density perturbation in the late 
Universe. 

As we have discussed above, e m is finite at r = and from Section III and IV, there is 
no long wavelength contribution to ( in the collapsing phase. It follows that e 2 = |eow^ 2 - ) in 
that phase. One situation of sepcial interest is where w ^ = in the collapsing phase, where 
the potential is irrelevant at <fi — > — oo and there is no radiation in the incoming state. In 
this case, £2(0^) = 0. In this case, we would obtain the same final result from any matching 
rule which set 62(0+) = v4e2(0~), with any constant A. 

The key point is that generically jumps across r = 0, since the background equation 
of state changes at the brane collision. Matching eo and €2 we obtain a long wavelength 
contribution to ( in the expanding phase, 

( lw (r > 0) ~ - A*;-2 eo ( w (2)< _ W W>\ (47) 
16 
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where -u/ 2 ^ and un 2 ^ are the values of for r < and r > respectively. Since e m oc A; 2 $ 
and, as discussed earlier, $ has a nearly scale invariant power spectrum, it follows from 
Eq. (|4~7|) that if undergoes a jump, then £ will inherit a scale invariant long wavelength 
piece. This is our main result. In the following sections we will study a simple example of 
a situation where is discontinuous. 

We shall need explicit formulae for the jump in w^ 2 \ and for €q. The formulae for 
are obtained from Eqs. and ( [Tip above. If we assume that prior to r = there is no 
radiation, so that in addition to scalar kinetic energy we have only the potential U, then 
taking the limit as r tends to zero from below one finds that 

(2)< _ Wl + IM^QU, tt ttIi t dU , dU 



u > = u ^ = a '^ +a '^ (48) 



where Q = ). Likewise, for the expanding phase, if we assume U = but 

radiation is now present, we obtain 

(2)> _ _8_(PrOT {m 

W 27 MW 1 J 

Since in general w^ 2 ^ 7^ w/ 2 ^, we infer that generically, a scale invariant spectrum of 
perturbations will, with our prescription above, propagate across r = into the expanding 
hot big bang phase. 

Finally, to compute the perturbation amplitude given in Eq. (|47| ) we need eo- This can 
be read off from the expression Eq. ( |l5l) above, by translating the r dependence into a'/a 3 
and employing the fact that the latter gives the exact dependence for the long wavelength 
modes of interest even when the potential breaks away from the pure exponential form used 
in the first half of this paper. We find that at r = 0, 

6 ° 3Mpi^/pa (a' Q - a^)' ^ 

To summarize the results of this section, we have elaborated the conditions under which 
a scale invariant spectrum of perturbations survives the passage through a = 0. Basically 
this requires that the equation of state, as parametrized by w(r), have a discontinuous first 
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or second derivative with respect to r at r = 0. This condition would appear to be quite 
generically fulfilled, in any situation where entropy is generated at a brane collision. A 
key assumption in the calculation was that the energy density is dominated by scalar field 
kinetic energy as we approach a = 0. This assumption is natural in the ekpyrotic scenario, 
where the scalar field represents the separation of the boundary branes. We shall explore 
one particular, simplified model in the next section. 

VI. BACKGROUND EVOLUTION IN A NEARLY LIGHT-LIKE BOUNCE 

In this section we review and extend the description of reversal from contraction to 
expansion, as elaborated in Reference 4. As discussed there, the background evolution is 
described by the variables ao and a\, which for the simplest brane model (i.e. branes in AdS) 
represent the scale factors on the positive and negative tension boundary branes. The four 
dimensional effective scale factor a is given by ■sJ and this vanishes at the bounce. 

If the potential V(<p) vanishes as runs off to — oo (more precisely, if the quantity a 4 V((j>) 
vanishes), then the Friedmann constraint equation implies that trajectory in the (ao, ai)- 
plane intersects the boundary of moduli space ao = a\ along a light-like direction!. Then, 
if no radiation is produced on the branes, the trajectory simply reverses, corresponding to 
the matching condition a' 01 (out) = — a 01 (in). We describe this as an elastic collision, since 
the internal states of the two branes are unchanged by the collision. 

However, at any finite velocity, the boundary brane collision must result in the production 
of radiation on the branes, since it is a non-adiabatic process. In the M theory context the 
corresponding string theory is weakly coupled near the collision, and this production of 
radiation should be computable once the correct matching conditions are understood. 

Let us consider the case where the incoming state has no radiation, and the potential V(<p) 
vanishes at t = and thereafter. This requires that the potential is turned off at collision, 
requiring a sudden and permanent change in the internal state of the branes. Associated 
with this change, we assume that a small amount of radiation, with density p r is generated, 
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on a brane with scale factor a. The collision is therefore inelastic. We parameterize the 
inelasticity as follows. The Friedmann constraint after collision reads 

a o(out) 2 - a; (out) 2 = A< f^ ■ ( 51 ) 

Since right hand side is positive, the outgoing trajectory must be time-like in the (a , air- 
plane. Since radiation redshifts as a" 4 , the expression (p^a 4 ) is a constant. If radiation is 
generated on both branes, this term can be taken to represent the sum of the corresponding 
terms for both branes. 

We then define the the efficiency £ with which radiation is produced by 

^ A(p r a A ) 

3M|,oi (in) 2 ' 1 ' 

which with fl5T|) yields a single equation for the two velocities a' (out) and a'^out). We need 
another equation to fix both. 

In the special case of dimensional reduction from five to four dimensions, there is a natural 
candidate for an approximately conserved quantity, analogous to the total momentum for 
an inelastic particle collision. As mentioned above, at small brane separations one has 
the standard Kaluza-Klein result that the size of the extra dimension is proportional to 
e y/V3^>/ M pi _ \ s stated above, we are assuming that the potential V(4>) vanishes as tends 
to — oo. Other terms in the Lagrangian describing matter on the branes may in principle 
acquire ^-dependence upon dimensional reduction. However, the terms describing massless 
gauge fields and fermions do not obtain any such ^-dependence due to their conformal 
invariance in four dimensions. (To see this, note that the four dimensional Einstein-frame 
metric g^ v is conformally related to the four dimensional components of the five dimensional 
metric.) Therefore at the classical level, in the <fi — > — oo limit, the Lagrangian describing 
gravity, <fi and four dimensional radiation possesses a global symmetry <p — > <p + constant, 
and it is plausible that the corresponding Noether charge Q, 

Q = -(cto^i — dia'o) oc a 2 (f)' , (53) 
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tends to a constant as the collision approaches. In this limit the classical equation of motion 
of 4> is just Q' = 0. 

However, the sign of Q must flip at the bounce. As argued in Ref. 4, this is essential in 
order that the trajectory remains in the physical region of the (ao, cii)-plane. The reversal 
of Q may be also understood by the following higher-dimensional argument. The value of 
Q at collision is proportional to the time derivative of log(ao/di). The latter quantity is 
the distance, rather than the vector displacement, between the boundary branes. Hence, if 
the branes come together and then draw apart, Q must change sign although its magnitude 
remains constant. It is therefore natural to impose the Z2 symmetry at t = 0, Q — > —Q, 
which should become exact in the limit that the collision velocity approaches zero. 

These arguments suggest that we parameterize Q- violation at the bounce (brane collision) 
using 

a _ Q(out) + Q(in) 

A = 7u~\ ' (54) 

Q(m) 



where A is expected to be small. Equations ([51]), (0) and (|5l ) together uniquely parame- 
terize the final values of a' and a[ after the bounce: 

od(out) = (l - A + 4 ^ - A) ) a ' l( ' in ' ) 

a', (out) = - - A - 4(1 1 A) ) a'i(in)- (55) 

An example of an ekpyrotic two-brane collision is shown in Figure 1, for a specific choice of 
the inter-brane potential. Both branes initially expand under the influence of the attractive 
potential. But when they get close, and the potential rises to zero, this decelerates ao so 
that it begins to contract. At collision, aQ (in) = — a'^in). Immediately afterwards, a' Q is 
positive and, for small £, a[ is negative. 
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FIGURES 




t 

(a) (b) 

FIG. 1. Sketch of a nearly light-like collision between two boundary branes. The potential 

employed was a 4 V((ft) = — af (oi/oo) 4 e - ^ where / = jq ((oq/oi) — 1)~ , chosen so that when 

expressed in terms of (ft, the potential V vanishes as eft — > — oo, in a manner mimicking the vanishing 

_ i _i 

of a non-perturbative potential 7 oc e ? or 7 oc e ». At collision, a and a[ are equal and 
opposite. The matching rule we propose is given in equation (p5|); in the Figure it is assumed 
that the efficiency £ and Q-violation parameter A are small. The scale factor oq is decreasing after 
collision: we assume it couples to a massless modulus \ which causes a\ to be repelled from a\ = 0. 
In the final state both ao and a± are expanding, and the universe becomes radiation-dominated 
with the outer-brane separation tending to a finite constant. 

Assuming that the potential remains zero after collision, and that only radiation is 
present, we have ao and a x flying apart linearly in t after collision, corresponding to an 
expanding a. According to (0), oq is decreasing after collision, and it tends to zero. This 
would lead to the separation between the branes going to infinity. However, it is easily 
avoided, essentially because the a\ modulus has a positive kinetic term. If there is one or 
more massless moduli fields x coupling to a\ as a\x' 2 '', they produce an effective potential 
in the a\ equation which is proportional to a{ 2 . This repels the a x modulus from cq = 0. 
Figure 1 shows an example of the full evolution, including the brane collision, reversal and 
turn-around of cq so that cto and cq are both expanding at late times. The final evolution of 
cq is insensitive to the value of x' i which only affects the evolution at small cq. In the long 
time limit, the scalar field (ft also tends to a constant and therefore so does the inter-brane 



28 



separation. 

One can generalize these considerations to examples where the potential U = a 4 V(4>) is 
finite and negative at brane collision. This requires rather special potentials V (</>), which 
diverge at large <f), but not too strongly. In this case, the trajectories are space-like at 
collision. 

Having specified the background evolution for the 'nearly light-like bounce', we can now 
consider the matching of perturbations. We have assumed that U vanishes at collision. For 
simplicity we shall assume that its first derivatives dU/dao and dll/dai also vanish there. 
Now we can read off from Eq. (|48 ) that u^ 2 ^ = 0, but we have from Eq. fl49|) and the 



c lw - C T - (58) 



formulae of this section that 

w(2)> = 2 aUinf _ 
3 a{ 

We can also read off from Eq. ( |50| ) that for the 'nearly light-like bounce', 

60 " ~ ZMp^a.a', (in)' (57) 
where we assumed £, A << 1. Putting these together in (H^), we find the final density 
perturbation amplitude is 

1 2 kr l ~ u a[(in ) 
12* M Pl ^p a\ 

It remains to compute a\ and a[ (in) at collision using the detailed behavior of the potential 
V{4>) for negative 0. 

As discussed in Ref. 4, as — > — oo, the string coupling constant tends to zero. It is 
natural to expect that the potential V(4>) goes to zero in this limit. We shall adopt a very 
simplified model here, in which the potential jumps to zero at some particular, negative 
value of <f)j. In this case it is straightforward to compute ai and a[ as r tends to zero from 
below. First, assume the jump happens at some time tj, before which equations (0) are 
valid. The total energy in the scalar field is + V(<j>) = 3p 2 M Pl /t 2 , and this is equal to 
the kinetic energy in the scalar field after the potential jump. So, just after the jump we 
have 
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Now, from a = 2acosh(0/^/6Mp ; ), we have 

a' = a (Hao - ^axj (60) 

and similarly for a^. Since the potential vanishes, and no radiation is present, a' and a[ are 
both constant up to collision. Then equations @, (|7|) and ( |6Qj) imply that 

a'o(in) = -2p(-t i ) 2p - 1 e^ /(v/5Afri) = -oi(in). (61) 

After the potential jump, ao and a% evolve linearly in r. To leading order in £ and A we 
have 

ao = (-tj) P (2cosh(^/v / 6MpO + ai(in)(r - r,)) 

oi = Hj) P (2sinh(-0 J / v / 6M„) - a 1 (in)(r - r,)) . (62) 
Setting these equal determines the time of collision and brane scale factors at collision, 

a = ai = (-t^Pe-fc/^w). (63) 

Now we have all we need to determine the final fluctuation spectrum. From equations (|58D, 
and fl63|), we have 



The dependence on k, tj and p is the same as that one obtains from the naive 'time delay' 
formula mentioned in Section II. Using equation (|16|) for 5<fi, the time delay method yields 
a perturbation amplitude ~ H5<j)/4> ~ Mp t \fpk~^ l+v ' {— t)~( 1+p \ in agreement with the 
dependence upon these quantities in equation (0). In the time delay argument, however, 
one uses the Hubble constant on the branes at collision, which is close to but not quite the 
same as the factor occurring in (|64|). 

Let us now translate the dependence of the last factor in equation into quantities 



determined by observations in the final expanding, radiation dominated Universe. After 
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collision, we assume the potential V{4>) is zero (because the internal state of the branes has 
changed so they no longer attract each other), so that the <fi modulus describing the inter- 
brane separation is a free massless field. We also assume radiation is present, in an abundance 
parameterized by £. It is straightforward to analytically solve the equations of motion post- 
collision, assuming the presence of the additional modulus needed to keep a\ away from zero. 
This modulus becomes irrelevant at late times. In the long time limit, one finds ao and a\ 
increasing linearly in conformal time r, with ao/ai — > |a' (out)/a / 1 (out)| ~ 1 + neglecting 
the dependence on A which is reasonable if A is small. Equating this to coth(— 0/\/6Mp/), 
we find that the final resting value for is given by 

However, what enters equation (|64j) is not <ftf, but <f)j, the value of <fi at which the potential 
V((ft) switches off. We can translate both values of <fi into the corresponding string cou- 
pling constants, using the relation g s oc eV^^ Mp! , which follows from M theory with the 
assumption that the six Calabi-Yau dimensions are fixed in the 11 dimensional metrical. 
For p close to zero, we find the final result 



^ ~\%%U(7)J • (66) 

The right hand side is the amplitude of the growing mode adiabatic density perturbation 
relevant to structure formation in the late Universe. (For an accurate calculation one should 
of course retain the p dependence, since over the many orders of magnitude of k involved, 
this can significantly affect the final normalization. We leave this complication for future 
work.) 

Our result for (j?" depends on the square root of the potential energy V at its minimum, 
in Planck units. This is reminiscent of the usual inflationary result. However, additional 
suppression factors arise. First, the numerical coefficient is small. Second, the factor £ 4 
is small if the efficiency of production of radiation at collision is small. Finally, the string 
coupling constant where the potential turns off, which we have crudely parameterized as 
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g s (j), would be expected to be substantially smaller than the value of the string coupling 
constant in the asymptotic outgoing state. Translating the formulae relevant to Horava- 
Witten theory, given by Witten,0 one finds for today's value of the string coupling constant 



where the volume of the Calabi-Yau manifold is VqMq UT . 

Before discussing numerical values, it is important to make the following caveats. First, 
the resting value of the scalar field <fi determined during the radiation era following the bounce 
is not necessarily that measured in today's universe. If there is a stabilizing potential for 
<f), that will instead determine the final resting value. Nevertheless it is conceivable that 
the resting value early in the hot big bang phase is closely related to the final value (as for 
example if develops a potential with many closely spaced degenerate minima). Second, 
the presence of additional moduli (such as are found in Horava-Witten theory) could have 
important consequences on the dynamical evolution of 4> in these early stages. In the above 
calculation we have limited ourselves to only one modulus, translating that directly into the 
string coupling constant. So the final numerical result can only be suggestive. 

For example, plausible values of the GUT coupling are gcur ~ 0.5, and the GUT mass 



10 17 GeV, giving g s {f) ~ 10 3 t> 2 . The turn-off of non-perturbative potentials might plausibly 
occur at g s (j) ~ 10, if instanton effects produce factors of the form exp(— 8tt 2 /g 2 ). The last 



factor in ( 16151) then yields ~ 10 _5 -u 3 . For v ~ 10~ 3 , and £ 4 (— Vj)/Mp t ~ 10~ 2 , we can 
obtain an amplitude ~ 10~ 5 , as required by observations. 

We conclude that the ekpyrotic scenario may offer a natural explanation for the smallness 
of the observed density perturbations. As we have emphasized, this is only suggestive at 
this stage, and will remain so in the absence of: 

• a microscopic check of the matching condition used for e m , within the context of M 
theory and string theory, 

• a computation of the efficiency parameter £ describing the production of radiation on 
the branes, 




(67) 
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• a check that the parameter A is indeed small, as was assumed, 

• and, a full calculation of a realistic inter-brane potential V(<j)) and a numerical solution 
of the equations improving the jump approximation used above. 

Although we have focused here on matching scalar perturbations, in principle one also 
needs a matching condition for tensor and vector modes. Neither acquire long wavelength 
power in the contracting phase of the ekpyrotic model, and it seems unlikely they will be 
generated at the brane collision. Nevertheless one can attempt to study possible matching 
conditions. It is not hard to see that the tensor amplitude hjj exhibits the same logarithmic 
divergence as the perturbation in the three-curvature of comoving slices, (. However, the 
canonically conjugate momentum a 2 hj- does tend to a finite constant at r = 0, suggesting 
it provides a possible matching variable. Again, establishing this will probably require a 
satisfactory microscopic theory. 

VII. CONCLUSIONS 

We have shown that inclusion of gravitational backreaction has a negligible effect on the 
density perturbations produced in the ekpyrotic universe, during the initial phase which is 
slowly contracting from the point-of-view of the four dimensional effective theory. 

We have proposed what we believe is a physically sensible matching condition at the 
bounce, based upon identifying physical perturbation variables which are well behaved (i.e. 
small) at r = 0, and matching on surfaces defined by the scalar field 0. The variables we 
used are well-behaved, we should emphasize, provided certain conditions are met as the 4d 
effective scale factor approaches zero, namely (1) the energy density is dominated by the 
kinetic energy of a scalar field (the modulus in our case) and (2) the inter-brane potential 
V{4>) does not diverge as tends to — oo. We showed that with the simplest matching 
prescription, namely matching the two linearly independent solutions at r = 0, the scale 
invariant spectrum of perturbations developed in $ early on in the contracting phase is 
generically passed on to the variable ( representing the amplitude of the long wavelength 
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growing mode density perturbation in the expanding phase. We also identified examples 
where no density perturbations are generated in the final Universe. If no radiation is gen- 
erated at the brane collision, and if the potential vanishes sufficiently smoothly there, the 
perturbations 'time-reverse' at a = and the amplitude of the growing mode perturbation 
is precisely zero in the final expanding Universe. Such examples are not realistic since there 
is no entropy generation at the brane collision, and in any case seem highly unlikely given 
that the outer-brane collision is not an adiabatic event. 

The existence of a 'zero-perturbation limit' is an intriguing feature of the ekpyrotic 
model and the matching prescription given here, since it suggests a natural explanation for 
the smallness of the observed density perturbations. Recall that a feature of inflation is that 
it naturally predicts a value of the density perturbation amplitude that is far too large, and 
fine-tuning of potentials is required to obtain a sufficiently small amplitude. 

We conclude that, while many issues connected to the microphysics at brane collision 
remain to be settled by rigorous investigation of string theory in the limit of outer-brane 
collision, the basic idea of Ref. [I] for producing density perturbations during the early stages 
of the ekpyrotic Universe, in a phase which is slowly contracting from the four dimensional 
point of view, remains viable. 
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Note Added: As mentioned in the Introduction, our conclusions differ from those 
of Lyth, Brandenberger and Finelli, and Hwang.i~i The disagreement is due to different 
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assumptions about physical conditions near the bounce, and to a different prescription for 
matching across it. We believe that the assumptions made by the authors of Refs. 7-9 are 
inconsistent with what we proposed in Refs. 1 and 4, because 

• they use a scalar field potential which diverges to minus infinity as the bounce is 
approached. This potential is not compatible with the bounce prescription discussed in Ref. 
4, since it is too singular. 

• they choose to match on a surface of constant energy density in the contracting phase, 
assuming that the scale factor of the universe reverses from contraction to expansion on this 
surface. This behavior is incompatible with the field equations describing the cosmological 
background solution. We instead follow the classical field equations all the way to the 
bounce, and apply a matching prescription consistent with our treatment of the background 
at this point. As we have explained above the appropriate matching surface in the ekpyrotic 
setup is defined by the scalar field specifying the inter-brane separation, and not by the 
energy density. 

• Brandenberger and Finelli claim our matching prescription based on the energy den- 
sity perturbation on slices of constant scalar field is inconsistent with results given in the 
literature for models where the equation of state undergoes a sudden jump0. It is easy to 
see why matching e m and its time derivative in this situation is incorrect. The equation 
of motion ( |33"1) for e m involves c 2 , related by ( |3~6l) to the time derivative of w. Hence, if 
w jumps, c 2 acquires a delta function contribution, which causes a jump in e' m across the 
matching point, which may be straightforwardly computed. The key point, however, is in 
the situation we are discussing, w is continuous across the bounce, thus there are no such 
delta function contributions. Therefore this 'counterexample' is truly a red herring. 

• Brandenberger and Finelli mistakenly imply that reversal from contraction to expansion 
would occur at a bulk brane-boundary brane collision. If the four dimensional effective 
description is valid, then reversal can only happen at a boundary brane-boundary brane 
collision, as described in Ref. 4. In the four dimensional effective description the universe 
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continues to contract after a bulk-boundary collision, so ( remains small and $ continues 
to grow. Only when the outer boundary brane collision occurs as it must, can the growing 
perturbations in $ get converted to long wavelength fluctuations in (. 
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